Nonfactorizable B — > Xco^ decay and QCD factorization 
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We study the unexpectedly large rate for the factorization-forbidden decay B — > XcoK within the 
QCD factorization approach. We use a non-zero gluon mass to regularize the infrared divergences 
in vertex corrections. The end-point singularities arising from spectator corrections are regularized 
and carefully estimated by the off-shellness of quarks. We find that the contributions arising from 
the vertex and leading-twist spectator corrections are numerically small, and the twist-3 spectator 
contribution with chiral enhancement and linear end-point singularity becomes dominant. With 
reasonable choices for the parameters, the branching ratio for B —> XcoK decay is estimated to be 
in the range (2 — 4) x 10~ 4 , which is compatible with the Belle and BaBar data. 

PACS numbers: 12.38.Bx; 13.25.Hw; 14.40.Gx 



B meson exclusive decays to hadrons with charmonium 
are interesting in studies of both strong interaction dy- 
namics and CP violation. The naively factorizable de- 
cays such as B -► J/tpK 0, 0, B -> r] c K 0, and 
B — > XciK H were studied in the QCD factorization ap- 
proach ||j in which the nonfactorizable vertex and spec- 
tator corrections were also estimated. 

To further explore the nonfactorizable contributions it 
is worth studying the factorization-forbidden decays such 
as B — > x c qK. Recently, B — > XcftK decay has been 
observed by Belle 0, and BaBar [3 with surprisingly 
large branching ratio which is even comparable to that 
of the factorization-allowed decay B — ► X&K : 



Bt(B + -> XcoK+) = (6.0±?:| ± 1.1) x 10- 4 [7], 

= (1.96 ±0.35 ±0.33) x 10~ 4 
Br(B ± -> XdoK*) = (2.7 ±0.7) x 10~ 4 [9]. 
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To explain the large decay rate of B — > XcoK, the final 
state re-scattering mechanism was suggested llOl. On 
the other hand, with the Light-Cone Sum Rules [ll| the 
nonfactorizable contributions were found to be too small 
to accommodate the observed B — > XcoK decay rate. 

In fact, in Ref. |a,ll2j within the QCD factorization ap- 
proach it was found that for the B — > XcoK decay, there 
exist both the infrared (IR) divergences in the vertex 
corrections and the end-point singularities in the leading 
twist spectator corrections. This implies that large non- 
factorizable contributions may come from soft gluon ex- 
changes. As argued in [12|, unlike the inclusive B decays 
to charmonium, where the IR divergences can be factor- 
ized into the color-octet matrix elements associated with 
the higher Fock states of color-octet cc with soft gluons 
[13 in nonrelativistic QCD (NRQCD) Q, the IR diver- 
gences in the exclusive two-body decays are difficult to 
be factorized. 

At the qualitative level, the results of Ref. ^3] may 
suggest that some fraction of the large color-octet con- 



tribution in the inclusive B decays to charmonium does 
in fact end up in two-body decay modes. That is, the 
soft gluon emitted by the color-octet cc pair may be re- 
absorbed by the quarks in the kaon or B meson, leading 
to a possible connection between the color-octet contri- 
bution and the infrared behavior of vertex and spectator 
corrections in QCD factorization approach. So qualita- 
tive estimates of these soft gluon contributions are im- 
portant for understanding the large branching ratios of 
both B — > XcoK and B — > J /^K . Furthermore, since 
the s quark emitted from the weak vertex moves fast in 
the B meson rest frame, we may expect that the soft 
gluon exchange is dominated by that between the cc pair 
and the spectator quark. 

In order to estimate the soft gluon contributions in 
these exclusive decays, we may let the charm quark be 
off the mass shell or give the gluon a mass, and then 
use the binding energy or gluon mass to regulate the IR 
behavior. Introducing the binding energy |15| or momen- 
tum cutoff [l(j was an useful way in estimating the inclu- 
sive annihilation rates of h c and Xd before the NRQCD 
factorization theory is developed. Even after that, this 
approach is still widely used in quarkoniu m p henomenol- 
ogy for qualitative estimates (see, e.g., though it 

is not a rigorous theory. Recently, the binding energy 
regularization was suggested in |l8j . In the present pa- 
per, we will estimate the non-factorizable decay rate for 
B — > XcoK by using the gluon mass regularization for 
the IR divergences, which is equivalent to the binding 
energy regularization in physical nature. The main dif- 
ferences between our calculations and those in Ref. 0] 
are the treatments of the end-point singularities in spec- 
tator interactions, which play the most important role in 
numerical evaluations. 

We treat charmonium as a color-singlet nonrelativistic 
(NR) cc bound state. Let p be the total momentum of the 
charmonium and 2q be the relative momentum between c 
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and c quarks, then v 2 ~ Aq 2 fa? ~ 0.25 can be treated as a 
small expansion parameter |14| . For P-wave charmonium 
XcOi because the wave function at the origin 7?-i(0) = 
0, which corresponds to the zeroth order in q, we must 
expand the amplitude to first order in q. Thus we have 
(see, e.g., [3) 



L Z ,S Z 



x^* JM (9)Tr[O"(0)P ls >,0)+O(0)Pf s > ) 0)] ) (2) 

where 0(q) represents the rest of the decay matrix el- 
ement and can be further factorized as the product of 
B — > K form factors and a hard kernel or as the con- 
volution of a hard kernel with light-cone wave functions 
of B and K mesons, within the framework of QCD fac- 
torization approach. The spin-triplet projection operator 
Pis z (p, q) is constructed in terms of quark and antiquark 
spinors as 1 

F is z (p,q) = \ —^v(~-q,82)u(~+q,siXsi;s2\lS.} 
V m c ^— ' 2 2 



ff J M N 



and 



(3) 

(4) 
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In Eq. @ we take charmonium mass M ~ 2m c in NR 
limit. Here m c is the charm quark mass. 

The integral in Eq. © is proportional to the derivative 
of the P-wave wave function at the origin 



^%«m(?) = *e* a (L z H/5<(0), (6) 



and we will use the following polarization relation for XcO 



l z s. 



X)e*^e*^<lL,;15,|00) = -L(-ff^+^).(ri 



Being contrary to Xco, the K meson is described rel- 
ativistically by the light-cone distribution amplitudes 



1 To construct the spin-singlet projection operator Poo(p,q), one 
only needs to replace f by 75 in Eq.JSJ. 
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where y and y = l—y are the momentum fractions of the 
s and d quarks inside the K meson respectively, and the 
chirally enhanced mass scale fix /(m s (/i)+md(£t)) 
is comparable to mj, which ensures that the twist-3 spec- 
tator interactions are numerically large, though they are 
suppressed by l/m^ (see 20]). The twist- 2 LCDA 4>K(y) 
and the twist-3 LCDA <fi p K (y) and <p^(y) are symmetric 
under y <-> y in the SU(3>) symmetry limit. In practice, 
we choose the asymptotic forms for these LCDAs, 



Mw)=0k(J/)= 61/(1 -»). = 1 



(9) 



The effective Hamiltonian relevant for B — > Xcoif is 
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Weff=%fK6K*iCl01+C 2 02)- Vt 6 ^Z)CiOi)-(ll»'l 

i=3 



where is the Fermi constant, Cj are the Wilson coef- 
ficients, and V qiq2 are the CKM matrix elements. We do 
not include the effects of the electroweak penguin opera- 
tors since they are numerically small. Here the relevant 
operators Oi are given by 

01 = (s a bf))v-A • {cf3C a )v-A, 

2 = (s a b a )v-A ■ (cpcp)v-A, 
03,5 = {Sgbg)v-A ■ / ] (QgilpW^A, 

046 = {s a bfi)v-A ■ ^2^(q0 a )v T A, (H) 

where a, [3 are color indices and the sum over q runs over 
u, d, s, c and b. Here {q\q%)v±A = 9x7^(1 ± 7s)<72- 

According to |f| all non-factorizable corrections are due 
to the diagrams in Fig^ These corrections, with opera- 
tors Oi inserted, contribute to the amplitude 0(q) in Eq. 
PJl. where the external lines of charm and anti-charm 
quarks have been truncated. Taking non-factorizable cor- 
rections in Fig^into account, the decay amplitude for 
B -^XcoK in QCD factorization is 



iM 



Gf 

7f 



VcbV^d-VtbVZid + Ce) 



x A, (12) 



and the coefficient A is given by 



A = 



^6^(0) a s C F 



(l- z )F 1 (M 2 )^-{fI + .fIl).{Yi) 



Here N c is the number of colors, C F = (N 2 - 1)/(27V C ), 
and Fq 1 are the B — > if form factors. We have used 



3 



where the fuctions d and c are written by 
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FIG. 1: Feynman diagrams for vertex and spectator correc- 
tions to B -^XcoK. 



the relation F {M 2 ) / Fi{M 2 ) = 1- z 0, Q, where z = 
M 2 /m 2 B « 4m;? /m^, to simplify the amplitude in (|13fl . 
The function // is calculated from the four vertex dia- 
grams (a, b, c, d) in Fig^ and fll is calculated from 
the two spectator diagrams (e, f) in Fig^ The function 
fll receives contributions from both twist-2 and twist-3 
LCDAs of the K-meson, and we can simply symbolize 
them as fll 2 and fll 3 , respectively. 

In the statement of color transparency 6] , the IR diver- 
gences should be cancelled between diagrams (a,b), (c,d) 
and (e,f) respectively in Fig^ But it is not true when 
the emitted meson is a P-wave charmonium, say, XcO- 
The soft gluons, which are emitted from a certain source, 
couple to the charm and anti-charm quarks through both 
color charge and color dipole interactions. For the color- 
singlet charm quark pair, the total charge is zero, but the 
color dipole interactions, which are proportional to the 
relative momentum q, give the leading order contribu- 
tions in both 1/mb power expansion and NR expansion 
(see Eq. @). As a result, there exist IR divergences in 
//, while fll 2, 3 suffer from logarithmic and linear end- 
point singularities if the asymptotic form of kaon LCDAs 
are used. 

Now we can express the function // as the Feynman 
parameter integrals 



fl=- 



d.r 



dy (d + c), 
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+ 



^xy((2z-6)x 2 + (6- zy-3y)x+y(z+2- zy)) 



((z + f)(*+f)) S 
y{3x - 2 + z{3y - 3x + 2)) 



+ 



z-1 
2 



K^-i^ + f) 

xy 2 (-3x + 2 + z(3x + y-3)) 



mz-l)x+f)¥ 



From Eq. 115(1 . we can see that the IR poles are all 
included in the function d before we integrate it com- 
pletely. The first two and the last two terms in d come 
from diagrams (a,b) and (c,d) in FigQ] respectively. The 
divergent integrals in Eq. I|14|) can be regularized by non- 
zero gluon mass m g , as we have mentioned above. And 
the gluon mass pole in the divergent integral is given by 



dx 
Jo 

iz(l 



dy d — 

-z + lnz) m 2 

- In ( — j) + finite terms. (16) 



(1-2T) 







Note that the infrared divergence would vanish if z 
(i.e., if treating the charm quark as a light quark). 

To derive functions fll 2 ' 3 , we extract the light-cone 
projector of K meson in momentum space from Eq. (|SJ), 



M Q >') = l -^{fj 5 <p K (y) 



d 



K 



(y) 



dk 2v 6 



)} fl . (17) 



where k%m is the momentum of the anti-quark (quark) in 
K meson, and the derivative acts on the hard-scattering 
amplitudes in momentum space. 

Using the projector given in i|17f) and eliminating 
(j) P x ' (y) by 10, we get the explicit form of fll 2 and fll 3 : 



fll 2 =a d£ 



fn 3 =- 



»(0 



dy $lM[-2z+(l-z)v], (18) 
V 



1 



di 



(0 



1 



10 ( Jo 1 

where the factor a is defined as 

87T 2 f K f B 



dy—[3z-(l-z)y], (19) 



a 



N c (l- z) 2 m 2 B F x (M 2 ) 



(20) 
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Here, due to our definition of fll in i|13fl . the form factor 
F\ is present in the denominator. 

In Eqs. (|18l 1191) . £ is the momentum fraction of the 
spectator quark in the B meson and 



rK(n) 



2m 2 K /{m b (n){m s (ij) + m d (n))} (21) 



is of order one and can not be neglected. The integral 
over f is conventionally parameterized as 



o s 



(0 



m B 
AT 



(22) 



If we choose the asymptotic form of the LCDAs of kaon 
in Eq. @, we will find logarithmic (linear) singularities 
in fll 2 (fll 3 ). As we have mentioned before, these singu- 
larities came from the color dipole interactions between 
the soft gluons and the P-wave charm quark pair, just 
like what happened in the vertex corrections. 

Our result of fll 2 is consistent with the previous ones 
[lil flij l. However, our function fll 3 is different from 
that in Ref. because they used a twist-3 light-cone 
projector 

m% = iJ r{f^My)-m7 5 ^-r K (y)} a0 , (23) 

which could be derived from Eq. (|17f) by adopting an 
integration by parts on y and dropping the boundary 
terms. If one simply parameterizes the linear singularities 
in Eq. {^J as 



0,2 



Ah ' 



(24) 



where ~ 500 MeV is the infrared cutoff, we will find a 
large difference between our results and those in Ref. . 
Since the decay rate is very sensitive to the numerical 
value of fll 3 as one can see in the following, we should 
be very careful in regularizing these singularities. 

The key point here is that the boundary terms have 
the form |* and can not be dropped safely, when 

the asymptotic form of </>J-(y) = 6yy is inserted. So the 
integration by parts is not well-defined when one simply 
parameterizes the linear singularities as in Eq. i|24|) and 
in Ref. That is why our result of fll 3 is different 

from that in Ref. 0] . 

The above problem comes from the observation that 
the product of virtualities of the propagators in Fig.l 
(e,f) goes to zero faster than <p^{y) in the end point 
regions. However, if we regularize all these small- 
virtualities carefully (i.e., introducing small off-shellness 
or transverse momenta for quarks and gluons), as sug- 
gested in Ref. 0], we can see that the boundary terms 
are exactly zero, 



y n r K { y) 
(y + A)' 



1+2 10 



0,1,2 



(25) 



and then the integration by parts is well-defined. That 
is, the two projectors are equivalent if and only if one 
regularizes the linear singularities properly, e.g., 



,/2 



ydy 
y 3 



dy = 1 
(y + A) 2 A 
ydy_ J_ 
" 2A 



(y + A) 



2 =--l + 0(A) 



1 + 0(A) 



(26) 



where the relative off-shellness A should be of order of 
Ah/niB- In practice, one needs to be careful and note 
that the integral kernels 1/y 2 and y/y 3 give different 
contributions in the scheme described by Eq. (|26p. al- 
though superficially they are the same in the expression 
of Eq. i|19fl and in the parametrization of Eq. (|24|) . 

It is worth emphasizing that this scheme is more phys- 
ical than the one in Eq. (|24|l since the off-shellness or 
transverse momenta of quarks and gluons are naturally 
serve as infrared cutoffs when y — > 0. Furthermore, it is 
a proper scheme to realize the factorization for the elec- 
tromagnetic form factors of tt at twist-3 level |22J . 

To determine the value of A, we introduce the bind- 
ing energy b ~ M — 2m c > [18j for x c0 , and then the 
virtuality of the quark line in Fig. 1(e) or (f) goes to 



(Pc + yp') 2 - m 



2 mB[V+ l-zM 



(27) 



Here it is evident that the last term in the last parenthe- 
ses plays the role of A. In the APPENDIX we will show 
that the relative off-shellness determined by the gluon 
propagator is of the same sign and the same order as A. 
In practice, we use 



A 



1-zM 



(28) 



Following the scheme in Eq. I|26|) and using the asymp- 
totic forms of LCDAs in Eq. 10, we reexpress the func- 
tion fll 3 as 



_ a-r K m B f 3z _ 3zy + (1 - Az)y 

l-z A B J Q m (y + X) 2 
3zy 2 - 3zy 3 



(y + A) 3 



(y + A) 4 



(29) 



For simplicity, here we have used the same relative off- 
shellness A to regularize each factor y in denominators. 
If we set A = in Eq. (|29|l . this expression will fall back 
upon the form given in Eq. I|19fl . We also derive the 
function by using the projector in Eq. (|23|l . and get a 
different expression, 



fTT 3 _ a - r K m B f 1 3z zy+(\-z)y 2 

; l-z'A B y dy[ (y + X) 2 (y + X) 3 



(30) 



Similarly, up to a universal normalization factor, the 
function fll 3 will fall back on the same form given in 
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Ref. |l8j if A is set to be zero. Completing the integrals 
in Eq. 1|29[) and Eq. (|30[1 . we get the same result 



fII 3 - a rK mB ' 5Z 



1 



I _ . Vy( 2X + (l-z)ln\+-(3-7z)+0(\)}, (31) 

as it should be. Again, we should emphasize that the 
emergence of infrared divergences and end-point singu- 
larities in decay amplitudes partly destroys the factor- 
ization assumption, and the soft interaction mechanism 
may be dominant in this decay mode. Before the diver- 
gences are removed or absorbed in some other factoriza- 
tion schemes, Eq. fpTTfl can reasonably serve as a model 
dependent estimation for soft gluon effects contributing 
to twist-3 spectator interactions. Furthermore, because 
of the linear singularities in Eq. (|31|l , fll 3 is not power 
suppressed in 1/mg, rather it is chirally and kinemati- 
cally enhanced. 

To be consistent with fll 3 , we regularize fll 2 in the 
same scheme and get 

fII 2 =a- 1 ^-[l2zln\ + 21z + 0{\)}. (32) 

For numerical estimates, we use the following input pa- 
rameters with two values for the QCD scale fi = yTOfcA/i ~ 
1.45 Gev and /i = m{,~4.4 Gev: 

M = 3.42 GeV,m B = 5.28 GeV, Xb — 300 MeV, 
/ B = 216 MeV 23], / A - = 160 MeV, Fi{M 2 ) = 0.75 H, 
^i(0) = VCK075GeV 5/2 [23, Ci(fi) = 1.239(1.114), 
C 4 (m) = -0.046(-0.027), C 6 (» = -0.068(-0.033), 
a a {p) = 0.34(0.22), r K {fi) = 0.85(1.3), A = 0.087.(33) 

In Eq. (|33|l the /i-dependent quantities at fi = 1.45 Gev 
(/i = 4.4 Gev) are shown without (with) parentheses. The 
Wilson coefficients Cj are evaluated at leading order by 
renormalization group analysis |2lj , since the amplitudes 
in Eq. i|13|) are of leading order in a s ■ 

The numerical results of fl are listed in Tab[l]with the 
gluon mass varying from 200 Mev to 500 Mev (the typical 
scale of Aqcd)- Comparing // with the values of fll 2 ' 3 , 
we find that for B — > XcoK, both // and fll 2 are small 
and they are partially canceled. As a consequence, the 
prediction for Br(_B — > XcoK) would be about an order of 
magnitude smaller than Eq. Q if only the leading-twist 
functions were used (see Tab[nJ. However, the chirally 
enhanced twist-3 contribution is numerically large and 
makes the predicted decay rate to be comparable to the 
experimental data Eq. Q. 

Our predictions for the branching ratio of B — > XcoK 
are listed in TablTTl and the values in the parentheses 
are the results evaluated by using leading-twist contribu- 
tions only. From TablTTl we can see that the results are 
not very sensitive to the value of the gluon mass. For 
comparison, we follow the available results in Ref. [Isj ]. 
where the IR divergences are regularized by binding en- 
ergy, to evaluate fl and give Br(_B^x cO i ; r)=22(16)xl0~ 5 



TABLE I: Functions evaluated by using the parameters in 
Eq. 1331 . The ^-dependent values given at fi = 1.45 Gev 
(fi — AA Gev) are shown without (with) parentheses. 





fl 


fll 2 


fii 3 M 


m g = 0.5 Gev 


12.9-17.7i 


-7.8 


35.0(53.6) 


m g = 0.2 Gev 


20.7-19.3i 


-7.8 


35.0(53.6) 



TABLE II: Theoretical predictions for the branching ratio of 
B — >XcoK including twist-3 contributions. For comparison, 
results without twist-3 contributions are listed in the paren- 
theses. 



10 5 x Br 


m g — 0.5 Gev 


m g = 0.2 Gev 


^ = 1.45 Gev 
M = 4.40 Gev 


30 (5.4) 
19 (1.7) 


42 (8.6) 
24 (2.7) 



for fi = 1.45(4.4) Gev. Again, the conclusion is similar 
to the gluon mass scheme. We see that by treating the 
singularities in spectator interactions properly, we can 
obtain the large experimental rate of B—^XcqK without 
introducing any unknown imaginary part for spectator 
interactions [18j. 

In summary, We have studied the factorization- 
forbidden decay B — * XcoK within the framework of the 
QCD factorization approach. We use the gluon mass to 
regularize the infrared divergence in vertex corrections. 
The end-point singularities arising from spectator inter- 
actions are regularized and estimated carefully by the 
off-shellness of quarks in the small virtuality regions. We 
find that for this decay the contributions from the vertex 
and leading-twist spectator corrections are numerically 
small, and the twist-3 spectator contribution with chiral 
enhancement and linear end-point singularity becomes 
dominant. With reasonable choices for the parameters, 
the B^XcqK decay branching ratio is estimated to be in 
the range of (2—4) x 10 -4 , which is compatible with the 
Belle and BaBar data. We would like to point out that 
there are other B exclusive decays to charmonia such as 
B->Xc2K, B^h c K as well as 5^(3770)^, and it is 
worth studying those decays to see if the soft spectator 
interactions are also the dominant mechanisms |26j. 

Note. Since this result was reported in 
arXiv:hep-ph/0502240, the B — > XcqK decay has 
also been studied with fcT-factorization in the PQCD 
approach 27], in which the vertex corrections are 
ignored and the spectator corrections are found to give a 
large enough decay rate in comparison with experiment. 
Their result is consistent with ours in the sense that 
in both approaches the vertex corrections are found to 
be small, and the spectator corrections give dominant 
contributions to the decay rate. For the vertex correc- 
tions, the authors of Ref. j2]j take them into account 
through the variation of renormalization scale for the 
factorizable contributions, while we use the gluon mass 
(or the binding energy) to regularize the IR divergence 



6 



and qualitatively estimate the size of these corrections. 
For the spectator corrections, in our approach the 
dominant contribution comes from soft gluon exchange 
or quark off-shcllncss, which are related to the end point 
singularities when the transverse momentum in the kaon 
is ignored. It is therefore interesting to examine if the 
main contributions in Ref. |2?| are also from the regions 
with small virtualitics for quarks. 
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APPENDIX 

In this APPENDIX we will show that the relative off- 
shellness determined by the gluon propagator in Figl.(e) 
or (f) is of the same order as A in i|28|) . To be explicit, we 
use the transverse momentum of the spectator quark to 
regularize the gluon propagator just like that in PQCD 
approach [isj ]. 

Let I denote the momentum of the spectator quark in 
B meson, then the denominator in the gluon propagator 
is (/ — &2) 2 — — (1 — z)t;ym B if the transverse component 
l± and k± are ignored. Then the transverse momenta can 
be used to regularize the gluon propagator near the end- 
point region y — > 0, and the denominator in this non-zero 
transverse momentum scheme can be written as |23l 



(I 



(1 - z)£ym 2 B - (l ± - k ± ) 2 



-(1 - z)m B {y + 



(1 — z)^m\ 



)• (34) 



serves as 



Obviously, here the average value , 

the relative off-shellness A for the virtual gluon. Note 
that all the components of I should be of order A — tub — 
Wfc, it is then easy to determine that A ~ 0(A/ms), 
which is of the same order as A in l|28|) . Furthermore, 
both A and A are of positive values, hence fll^ 2 ' 3 ^ in 
Tab. I do not contain imaginary parts. 
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